Screening effects in Coulomb frustrated phase separation 
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We solve a model of phase separation among two competing phases frustrated by the long-range 
Coulomb interaction in two and three dimensions (2D / 3D) taking into account finite compressibility 
effects. In the limit of strong frustration in 2D, we recover the results of R. Jamei, S. Kivelson, 
and B. Spivak, Phys. Rev. Lett. 94, 056805 (2005) and the system always breaks into domains in a 
narrow range of densities, no matter how big is the frustration. For weak frustration in 2D and for 
arbitrary frustration in 3D the finite compressibility of the phases is shown to play a fundamental 
role. Our results clarify the different role of screening in 2D and 3D systems. We discuss the 
thermodynamic stability of the system near the transition to the phase separated state and the 
possibility to observe it in real systems. 

PACS numbers: 71.10.Hf, 64.75.+g, 71.10.Ca 



I. INTRODUCTION 

In the presence of long-range forces the tendency to- 
wards phase separation among two competing phases 
may lead to self-stabilized domain patterns of finite 



size 



1.2.3.4.5 



A prominent example is provided by charged 
systems compensated by a rigid background. The 
Coulomb energy cost grows faster than the size of 
the system and precludes macroscopic phase separa- 
tion. Therefore one finds either uniform phases or 
phases with domains of one phase hosted by the other 
one^&iS&i ^! 2 ^! In the mixed state the finite 
Coulomb energy arising from the mismatch between the 
local electronic charge — en^ and the fixed background 
charge en, competes with the surface energy to deter- 
mine the typical size of the domains. 

Experiments in layered materials like cuprate o 15 ! 16 ' 17 , 
ruthanatesi^, manganite thin films*^ and in the two- 
dimensional electron ga o 20 i 21 have fueled interest in the 
two-dimensional version of this problem— i 14 i 22 . In par- 
ticular the recent discovery of pronounced anisotropics 



m transport propertie; 



,18,20 



are consistent with the pro- 



posal of exotic electronic liquid phases 23 analogous to the 
intermediate order states of liquid crystals 24 . 

Jamei, Kivelson and Spivak 13 have analyzed a model in 
which the energies of two infinitely compressible uniform 
phases cross each other at a certain density n c . They 
have shown that in 2D, sufficiently close to n c , the sys- 
tem always breaks into domains for any level of frustra- 
tion due to the long-range Coulomb (LRC) interaction. 
They have focused on the universal aspects arising for 
large frustration where the compressibility turns out to 
be negligible. On the other hand it has been argued that 
screening, driven by finite compressibility, plays a funda- 
mental role on Coulomb frustrated phase separations^. 
It is thus of interest to analyze its effect for the full range 
of frustration and, in this respect, to clarify the difference 
between the two-dimensional (2D) and three-dimensional 



(3D) case. 

In this work we solve a model of Coulomb frustrated 
phase separation, similar to the one of Ref. [13| but aug- 
mented to take into account the finite compressibility of 
the phases (Sec. |TT|. This allows us to derive the phase 
diagram in 2D from the limit of strong frustration down 
to the limit of zero frustration where phase separation is 
determined by Maxwell construction (MC) and the com- 
pressibility of the phases, by sure, cannot be neglected. 
The interplay between the screening and the size of in- 
homogeneities is discussed in two-dimensional systems 
(Sec. IIII[) and three-dimensional systems (Sec. lIV|) . Both 
cases are compared at the end of Sec. HVl We show 
that for not so large frustration in 2D and practically 
everywhere in 3D the physics of the mixed state can be 
captured in a simple approximation where the density 
variation inside the inhomogeneities is neglected. We ar- 
gue that the strongly frustrated 2D mixed state regime 
is experimentally hardly accessible and characterized by 
a large negative electronic compressibility that can lead 
to a volume collapse of the lattice. The model presented 
here can be solved exactly and confirms the conclusions 
of previous approximated treatments. In particular it 
illustrates the general rule found in Refs. [gflll that for 
generic value of the parameters inhomogeneities can not 
have all linear dimensions larger that the 3D screening 
length. 



II. MODEL AND GENERAL SOLUTION 



We consider a ferromagnetic Ising model linearly cou- 
pled to a charged fluid as a generic model of a first order 
density driven phase transition with Coulomb interac- 
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FIG. 1: (Color online) Behavior of the free energy density for 
the uniform phases with s = ±1 in the present model (thick 
lines) compared with the ft — > oo model discussed in Reffl3l 
(thin full lines). The dashed lines are the energy density of 
the mixed state solutions in 2D systems for different values of 
the frustrating parameter defined in Sec. lIIII The energy den- 
sity is measured in unit of the characteristic phase separation 
energy gain en = (5n°) 2 /K = 2irQ 2 (Sn°) 2 l s . The dot-dashed 
line corresponds to the Maxwell construction (A = 0). 



tions. The Hamiltonian reads: 

H = —J ^ SjSj ^ 3j (nj - n c ) (1) 

<ij> i 




where Sj = ±1, the index i runs over the sites of a hyper- 
cubic lattice of dimension D — 2,3 and n is the average 
charge. Instead of the present Ising term one could con- 
sider, as in Ref. [HI a double- well potential in the very 
large barrier limit. In solid state systems the strength 
of the Coulomb interaction will be given by Q 2 = e 2 /eo 
with eo the static dielectric constant of the environment 
(i.e. excluding the mobile electrons). 

For the two uniform phases (rij = n), the energy as a 
function of density consists of two parabolas that cross 
at n c with chemical potential difference A/x c , minima at 
n = n c ± Sn° /2 and Sn° = kA/i c as illustrated by the 
thick lines in Fig. [1] For Q = 0, according to the ordi- 
nary Maxwell construction, the system phase separates 
when the average density ri is in a window of width Sn° 
around the crossing density. The energy of this solution 
is depicted with the dot-dashed horizontal line in Fig. Q] 

We are assuming that the two uniform phases have 
the same compressibility to make the problem solvable. 
A related model in 3D has been considered in Ref. 0. The 
2D case in which one of the phases is incompressible has 
been analyzed in Ref. [HI, whereas Ref. [lU corresponds to 
the case in which both phases are infinitely compressible 
and Maxwell construction is not defined (n — > oo, thin 
full lines in Fig. []J . 

Interfaces of the Ising order parameter are sharp 



by construction with a surface tension given by a — 
2J/a D - 1 where a is the lattice constant. The linear term 
in A/j, c expresses that the homogeneous phase with Si = 1 
(t) is favored at higher densities, and the phase with 
Si = — 1 (|) at lower densities. 

We solve the problem in the limit in which domains 
are much larger than the lattice constant and replace the 
site index by a continuum variable, i.e. s(r), n(r) in a 
volume V. Hereafter we take a = 1. For each mixed 
state, the geometry of the domains defines s(r) and its 
Fourier transform s(q). Writing the energy in the Fourier 
space, the q = component of the Coulomb term is 
canceled. At q ^ 0, upon minimizing with respect to the 
charge distribution, one gets: 

n(q) (k- 1 + 2 D -\Q 2 /\ q \ D -i) = ^s(q) (2) 

This equation has a simple physical interpretation by 
noticing that for a fixed domain configuration 

M M = y s W + I n W ~~ " c l 

represents the local chemical potential and the electro- 
static potential can be put as: 

v (r) = J dr'cj) (r — r') n (r') . 

Here <j) (r) is the Fourier transform of 

as follows from Poisson equation in three and two- 
dimensional systems^ in the presence of the 3D Coulomb 
interaction and plays the role of the "effective" Coulomb 
interaction. With these definitions Eq. ([2]) states that 
the total local electrochemical potential fj, (r) + v (r) is 
constant. The latter condition is the generalization to 
electronic systems of the Maxwell condition for neutral 
fluids that is enforced by the constancy of the local chem- 
ical potential across different phases. 

We can now use Eq. ([2|) to eliminate the charge from 
the energy and obtain an energy functional that depends 
only upon s(q): 

E = all ^s (n - n c ) + ^- (n - n c f + (3) 

Z Ik 

A^c \ ^ s q s -q 

sv ^ o [«- 1 + 2 c -%g 2 /| q p-i] 

where we have dropped an irrelevant constant. The first 
term in Eq. ([3]) is the surface energy of the sharp inter- 
faces with S indicating the total surface of the interfaces. 
The second and third terms are the q = contribution 
from the bulk free energy of the competing phases. The 
last term comes from the q ^ contribution of the last 
three terms of Eq. |T]) after eliminating the charge via 
Eq. ©. 
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FIG. 2: (Color online) The phase diagram for the smectic 
stripe solution in 2D. The central region corresponds to the 
mixed state (MS). The thick line is the exact solution of the 
present model with finite compressibility, the thin full line is 
the hi — ► oo limit taking the short distance cutoff of Ref. [l3| 
equal to the present screening length l s . The dashed line is 
the UDA. The inset shows an enlargement near the crossing 
density n c where the UDA predicts a critical value of the 
frustrating parameter. 



Eq. has to be optimized with respect to the ge- 
ometry of the domains. Here we restrict to a periodic 
structure of alternating s = ±1 stripes. This state 
can be interpreted as a smectic electronic liquid phase 
that possesses orientational order and breaks the trans- 
lational symmetry only in one direction. Possibly, fluc- 
tuations of the stripe order can restore the translational 
symmetry 2 ^ thus leading to a nematic phase. For the 
smectic stripe state the problem can be solved analyti- 
cally (Sees. UlIIIVp . Other states with different geome- 
tries like circular drops can be treated numerically and 
will not be considered in this work. 



III. 2D SYSTEMS 

In this section we consider 2D systems. From Eq. ([3]) 
the energy density (e = E/V) of a mixed state with 
alternating s = 1 stripes of width 2Rd and s = — 1 stripes 
of width 2(R C — Rd), can be cast as: 
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where we introduced the two-dimensional screening 
length 

l s = (2nQ 2 K)-\ 

and the volume fraction v = Rd/R c of the s = 1 phase. 
We also define the Fourier transform of the effective in- 
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with n running over nonzero integers and 



\q\(l + h\q\Y 

This effective interaction coincides in this case with 
the effective 3D screened Coulomb interaction in two- 
dimensional electronic systems discussed in Refs. and 
0. 

Eq. ((3]) can be recast in dimensionless form measuring 
the energy in units of the characteristic phase separation 
energy density gain eo = (Sn ) 2 / k, lengths in units of l s 
and densities in units of Sn°: 



e 
eo 



— = A 



2R C 



v 

2 



2R C 



n — n c 
5n° 

9 Rc 



Sn° 



(5) 



Here we defined a dimensionless coupling: 

A 2 = 4trtQ 2 / A^c 

that parametrizes the mixed state energy density. Apart 
from numerical constants, A coincides with the frustrat- 
ing parameter defined in Refs. IMlil that measures the 
strength of the frustration due to LRC interaction and 
surface energy with respect to eo- 

The dimensionless free energy density, Eq. ([5]), depends 
upon A, (n — n c )/5n°, v and R c jl s . Minimizing with 
respect to the volume fraction v and R c /l s we obtain the 
phase diagram in the A-density plane shown in Fig. [5] 
with the thick lines. 

As in many other situations, increasing the frustration, 
the region of stability of the uniform phases is increased. 
As for the n — oo case (thin full line in Fig. [2|), there 
is no direct first-order transition between the two homo- 
geneous phases no matter how big is A^. Indeed the 
transition line is logarithmically singular at n = n c and 
thus for large A the mixed state is enclosed in the expo- 
nentially narrow range 



\n — n c \ < 5n e 



(6) 



Clearly, observation of the mixed state for large A re- 
quires an enormously accurate control of the density n 
which may be hard to achieve in practice. As frustration 
decreases, the range of densities of the mixed state grows 
and converges to the Maxwell construction range when 
A^O. 

Close to the divergence, the present phase diagram co- 
incides with the one for infinite compressibility (thin full 
line in Fig. O provided one takes the short-distance cut- 
off of Ref. [l3| equal to l s . This agreement is due to the 
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FIG. 3: (Color online) Charge density modulation for a cut 
perpendicular to a stripe in 2D at n = n c , Rd = Rc/2 and 
three different values of the screening length corresponding 
to A ~ 0.2 (R c = l t ), A ~ 0.7 (Re = Wl s ) and A ~ 1.1 
(R c = 100/s). In the case R c /l s — > oo the charge density 
diverges at the boundary. A finite l 3 removes the divergence. 
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FIG. 5: (Color online) Behavior of the width (top panel) and 
periodicity (bottom panel) of \ stripes in 2D systems for dif- 
ferent values of A < 1. The dashed lines correspond to the 
uniform density approximation whereas the full line is the ex- 
act solution. Notice that the \ stripes are the minority phase 
for n < n c and the majority phase for n > n c 



FIG. 4: (Color online) Behavior of Rd (log scale) at the cross- 
ing density as a function of the frustrating parameter. The 
thick line is the solution of the present model while the thin 
line is the infinite compressibility limit* 3 . The UDA is the 
dashed line. 

fact that in the highly frustrated regime the physics is de- 
termined by the slow power-law relaxation of the charge 
(Fig. [3]) far from the stripe boundary. In fact the stripe 
width is much larger than the screening length and in- 
deed behaves as (c.f. Fig. 2]) 

R d = ^~l s e* x2 (7) 

whereas the finite compressibility of the present model 
affects the charge only in a range of order l s from the 
boundary. Its effect is to remove the unphysical diver- 
gence of the charge density at the stripe boundary arising 
when l s = (Fig. [3]). . On the contrary, for low A the 
stripe size is of the order of the screening length Rd ~ ^A 
(Fig. 3]) and the finite compressibility is relevant. 



For practical applications and for more complicated 
forms of the energy specifying the homogeneous phases, 
it may be convenient to use a uniform density approxima- 
tion (UDA) in which one neglects the spatial variation of 
the density inside each inhomogeneity^i^. The phase di- 
agram in this case is also plotted in Fig. [2] with a dashed 
line. As in Refill the electrostatic energy has been ap- 
proximated by taking into account only the self inter- 
action i.e. we neglect the Coulomb interaction among 
different stripes. 

Remarkably the UDA with this crude approximation 
for the Coulomb interaction gives very accurate results 
in a wide range of the phase diagram except around the 
crossing density, since it misses the logarithmic singular- 
ity and then predicts a critical value A c of the frustrating 
parameter above which the uniform phases would be sta- 
ble as in the 3D case^ and a first order phase transition 
among them would occur (see the inset of Fig. [2]). Ap- 
proaching A c from below, the MS disappears with a di- 
vergent stripe size (see Fig. H]) in contrast with the exact 
results where the MS persists with a finite stripe size. 

For A < A c , at the onset of the mixed state the UDA 
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FIG. 6: (Color online) Behavior of the period of layers (log 
scale) vs A at the crossing density in 3D (thick line) and 2D 
systems (dashed line). 



predicts a first-order transition with the minority phase 
stripe size and the stripe periodicity that stay finite (see 
Fig. [5]). Taking into account the local dependence of 
the electronic density, the jump in the periodicity of 
the stripes at the transition is substituted by a square 
root divergence. Thus the uniform-mixed state transition 
is characterized by a finite value of the minority phase 
stripe size and a divergent periodicity which corresponds 
to a second-order transition^. Apart from the singular- 
ity at the transition between the uniform phase and the 
smectic stripe phase, the behavior of the mixed state for 
not so large frustration is very well represented by the 
UDA (Fig. [HE]). 

For A > A c relaxation of the UDA produces the loga- 
rithmic singularity which, as stated above, affects in an 
exponentially narrow region the phase diagram. Further- 
more, according to Eq. ((7]), for large A the size of the do- 
mains can quickly reach the size L of the system. Finite 
L will introduce a cutoff to the singularity 

A:~[(logL/Z s )/7r] 1/2 . 

For A > A* the transition will look first order like. 

For A < A* the transition is second order like, however 
the effect of the mixed state solution is to produce an 
exponentially small rounding of the singularities in ther- 
modynamic quantities quite hard to distinguish from a 
first order jump. On the other hand an important differ- 
ence from a true first order transition is that the latter 
shows hysteresis when driven at a finite rate whereas here 
hysteresis will be absent for A < A* due to the fact that 
the surface energy is effectively negative. For A > A* 
instead, hysteresis can appear. 

Now we discuss the thermodynamic stability. In Fig. [1] 
we show the energy of the mixed state solutions. For fi- 
nite A the energy vs. density has a negative curvature 
in the mixed state indicating a negative electronic com- 
pressibility. For large A the electronic compressibility in 



^/l - [e^^n- n c )/6n' 



012 



(8) 



kJ 1 negatively diverges at the uniform- mixed state tran- 
sition and is exponentially large and negative for n ~ n c . 

A negative divergence of the compressibility with a 1/2 
critical exponent also arises at the uniform-stripe transi- 
tion for small A but with a strength that vanishes when 
A -> 0. 

In the present model the background is assumed to be 
incompressible (fc& = 0) but real systems will have a small 
background compressibility fcj, > and the system will 
become unstable when the total inverse compressibility 
k^ 1 + k^ 1 < 0. This will lead to a volume instability 9 , 
substituting the stripe transition with a volume collapse 
transition and reintroducing hysteresis 28 . This is the 
most likely behavior to be found in real systems specially 
for large A where the electronic compressibility is large 
and negative in the whole mixed state stability range. 



IV. 3D SYSTEMS 

In this section we discuss the 3D case. The energy of 
a layered mixed state is given by Eq. (U) with D = 3, 
provided the effective interaction u (q) is given by 

u(q) = (l + iy)-\ 

and the screening length l 8 is substituted by the usual 
3D expression, 

l s = (4ttQ 2 K )- 1 / 2 . 

Notice that the standard 3D screened Coulomb interac- 
tion is given by u (q) /q 2 . 

The different nature of the screening in two- and three- 
dimensional systems affects strongly the properties of the 
mixed state. At the crossing density, the exact expression 
for the energy density of a layered state is minimized for 
v = Rd/R c = 1/2 an d takes the following simple form: 



(7 



I-A- 3 / 2 tanh^ 



(9) 



where A is the 3D frustrating parameter; 9 
A = 4[(,QV)/( K A^)] 1/3 . 

The energy density of the uniform state is given by 
e = 0, hence for a layered state to be possible the term 
in the brackets in Eq. §§§ must become negative. This 
condition is equivalent to the existence of a critical frus- 
trating parameter A c = 1 as it was derived in Ref. 
within the UDA and now showed in an exactly solvable 
model. 
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FIG. 7: (Color online) Charge density modulation for a cut 
perpendicular to the layers in 3D systems (thick line) and 2D 
systems (dashed line) at n = n c and a layer period R c — 20l s . 



Minimizing Eq. §§§ respect to R c one finds that R c is 
smaller than a few screening lengths except in the un- 
physical case in which A is fine-tuned exponentially close 
to A c . This is in contrast with the 2D case where R c is 
unbounded for a generic large A (Fig. [S]). 

This result is due to the different screening effect in 
three-dimensions with respect to the 2D case. To illus- 
trate this difference we plot in Fig. [7J the charge density 
in 2D and 3D in a cut perpendicular to the stripe/layer 
for R c = 20^. In 3D the charge density decays exponen- 
tially from the interface over a distance of the order of the 
screening length whereas in 2D it decays as a power law. 
As emphasized in Ref. 9 the phase separation energy gain 
stems from the regions where the local density n(x) is sig- 
nificantly different from the global n value. In 2D this is 
generically fulfilled and thus phase separation is favored. 
In 3D instead, the two densities are substantially differ- 
ent only in a region of width l s around the interface. For 
A < A c the phase separation energy gain from these re- 
gions compensates the surface and Coulomb energy cost 
and makes inhomogeneities possible. Instead the central 
region in the 3D case of Fig. [JJ produces an exponen- 
tially small phase separation energy gain and therefore 
one never finds inhomogeneities with R c >> l s , as in the 
figure, unless A is exponentially close to A c . Indeed one 
can check from Fig.[S]that this is the case for the present 
value of R c /l 8 . 

The fact that strong variations of the density as the one 
pictured in Fig. [JJ practically never occur in 3D makes the 
UDA reliable everywhere in the 3D phase diagram. 

The exact 3D phase diagram for the layered state in 
the A-density plane is shown in Fig. [51 As expected it 
is very similar to the one obtained within the UDA in 
aflU The only difference is that in the UDA the 



Refs. 

uniform-inhomogeneous transition resulted weakly first 
order whereas here it is second-order like due to a diver- 
gence of R c /Rd as in 2D. 

The main difference with the 2D case is that the two 



FIG. 8: (Color online) The phase diagram in the A — (n — 
n c )/8n° plane for three-dimensional systems. The thick lines 
correspond to the second-order transitions from the homoge- 
neous phase to the smectic stripe phase characterized by a 
divergent periodicity and a finite minority phase layer size. 
The thin line represents the first order transition between the 
two homogeneous phases that occurs for A > A c . 



second order transition lines from the homogeneous phase 
to the mixed state touch each other at A = A c . For 
A > A c the systems has no intermediate state between the 
uniform phases and a direct first-order transition appears 
among them. In this case the energy density as a function 
of the global density is given by the lowest energy branch 
among the two parabolas shown in Fig. [TJ This has a 
cusp singularity at n c which produces a Dirac function 
like negative divergence of the compressibility. This leads 
to a volume instability also in this case&2£. 



V. CONCLUSIONS 

In this work we have solved a model of phase separa- 
tion frustrated by the LRC interaction while considering 
screening effects and the charge relaxation inside the do- 
mains in two and three dimensional systems. 

In 2D and for high frustration the finite bare com- 
pressibility can be neglected. Its effect can be taken up 
by a short distance cutoff. In this regime, however, the 
mixed state requires very stringent conditions to be ob- 
served: i) it requires an enormously accurate control of 
the density since it appears in an exponentially narrow 
region around the crossing density, ii) the size of the do- 
mains grows exponentially with the frustration and may 
easily reach the size of the system , Hi) the mixed state 
is characterized by an exponentially large negative com- 
pressibility which induces a volume instability when the 
rigidity of the background is not sufficient to stabilize the 
system altogether. 

An incompressible background may seem unphysical 
according to the above discussion. However in some cases 
a sufficient separation of energy scales may avoid a vol- 
ume instability if singularities in the compressibility are 
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rounded by extrinsic effects. For example in ruthanates 
and in the 2D electron gas the relevant electronic phe- 
nomena occurs at temperatures below a tenth of a Kelvin 
to be compared with melting temperatures of the ma- 
terial of the order of hundreds of Kelvin. Even when 
a volume instability may be avoided, the strongly frus- 
trated mixed state may appear somewhat academic due 
to the stringent conditions on the density. An important 
physical, however, consequence is that the transition will 
look first order like but without hysteresis as discussed 
in Sec. IIIII In the case of the ruthanates there is the 
additional advantage that the control parameter is the 
magnetic field which allows for considerable fine-tuning. 

For low frustration the finite bare compressibility of the 
phases [the third term in Eq. (p}] cannot be neglected. 
The mixed phase compressibility gets strongly renormal- 
ized. We have discussed its behavior in 2D [Eq. fSJ] 
which in principle can be directly measured by capac- 
itive techniques^. We have also shown that neglecting 
the spatial variation of the charge inside the domains, one 
obtains very accurate results at low frustration. From 
our results and Figs. 12141 it is clear that the limit k — > oo 
and the UDA can be seen as the strong and the weak 
frustration approximation respectively for the full model 
Eq. p. 

In 3D systems the bare compressibility term is essential 
to obtain the mixed state with sharp interfaces as studied 
here and the UDA becomes reliable everywhere in the 
phase diagram. 

The different role of the screening in 2D and 3D can be 
interpreted by noticing that the electronic charge density 
behavior inside the domains is dramatically different in 
3D with respect to 2D. In the latter case, in fact, the 
charge density inside one stripe decays with a power-law 
from the interface and the local electronic charge density 
differs from the background density (i.e. n (x) ^ n) over 
the entire stripe width. Thus domains gain phase sepa- 
ration energy in all the region and there is no limitation 
for their size. This also explains why the mixed state can 
appear independently on how strong the frustration is. 

On the contrary, in 3D systems the charge density 



decays exponentially. As a consequence, the system is 
forced to satisfy a "maximum size rule'^ill that states 
that for generic parameters the inhomogeneities can- 
not have all linear dimensions much larger than the 3D 
screening length . This allows for arbitrary large inho- 
mogeneities in 2D since one of the dimensions is already 
smaller than P S D as indeed found (Fig.0. 
The frustrating parameter can be written as 



>D-1 



Here Id represents the size that inhomogeneities should 
have for the surface energy cost, a /Id, be of the same or- 
der as the phase separation energy gain, eo = (dn°) / k. 

To make a rough estimate of A we go back to the lat- 
tice and take the bare compressibility as an inverse band- 
width k — [2Dt] 1 with t an interatomic hopping inte- 
gral. Defining a nearest neighbor Coulomb interaction 
V = Q 2 /a one obtains: 



A 2 
A 3 



7T VJ 



2t 2 (5n ) 2 
32tt VJ 2 
27 t 3 (<5n ) 4 



(D = 2) 
(D = 3) 



If one neglects the logarithmic singularity in 2D, inhomo- 
geneities require A < 1. In general we expect the scale 
J which may be due to magnetism or other low energy 
phenomena to satisfy J < t. According to our definitions 
Suq is the difference of the Maxwell construction densi- 
ties (measured as number of particles per site) . Typically 
we expect weak phase separation tendency so Suq << 1. 
Comparing with the 2D case, the 3D case has an extra 
factor J/t that makes the mixed state more favorable 
and an extra factor Sn^ 2 that makes the mixed state less 
favorable. A more precise analysis requires microscopic 
modeling of specific situations and will be presented else- 
where. 
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